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Many relatively thick metal oxide films grow according to what is called the parabolic law L — 
y/2At + . . . . Mott explained this for monovalent carriers by assuming that monovalent ions and 
electrons are the bulk charge carriers, and that their number fluxes vary as t~^^^ at sufficiently 
long t. In this theory no charge is present in the bulk, and surface charges were not discussed. 
However, it can be analyzed in terms of a discharging capacitor, with the oxide surfaces as the 
plates. The theory is inconsistent because the field decreases, corresponding to discharge, but there 
is no net current to cause discharge. The present work, which also includes non-monovalent carriers, 
systematically extends the theory and obtains the discharge current. Because the Planck-Nernst 
equations are nonlinear (although Gauss's Law and the continuity equations are linear) this leads to 
a systematic order-by-order expansion in powers of t~^^^ for the number currents, concentrations, 
and electric field during oxide growth. At higher order the bulk develops a non-zero charge density, 
with a corresponding non-uniform net current, and there are corrections to the electric field and the 
ion currents. The second order correction to ion current implies a logarithmic term in the thickness 
of the oxide layer: L = V2At + B Int + . . . . It would be of interest to verify this result with 
high-precision measurements. 



I. INTRODUCTION 

From the late 1930s to the late 1940s, N.F. Mott 
considerecPfS the implications of experimental resultJ^ 
for oxide growth. (For more of a review, see Ref [51) Un- 
der certain conditions (particularly high temperature), 
many metals develop - a layer of oxide at a parabolic rate 
on surfaces exposed to gas containing oxygen, according 
to 



= 2At, 



(1) 



where L is the thickness of the oxide, t is time, and A is 
a constant. The rate of growth is thus 



(2) 



This result may be thought of as representing the long- 
time asymptotic limit. For specificity, we assume that 
the metal M fills a; < 0, that oxide MO fills < x < L, 
and that oxygen gas O fills L < x. This means we employ 
a moving coordinate system where x = represents the 
M/MO interface, and x ~ L{t) represents the MO/0 
interface. 

A field and fluxes that vary as t^^/^ are expected on 
the basis of a gradient of concentration, with the val- 
ues of the carrier concentrations pinned by the two sur- 
faces and the length L determining the gradient.'' That 
is, dL/dt ^ 1/L gives a parabolic law. Wagner obtained 
a parabolic growth law using the Planck-Nernst equa- 
tions and some additional assumptions.'^ Mott obtained 
a parabolic growth law using a more complete argumenlP 
that invokes the Planck-Nernst equations. Gauss's Law, 
and (implicitly) the continuity equations. In this case one 
can think of the electrochemical potentials pinned by the 
two surfaces and L determining the gradient, which then 
yields the parabolic law. 



For electron and ion number currents (jaijb) and ion 
valence Z ^ 1, Mott assumed that the total current J ~ 
e{jb — ja) in the oxide is zero, so 



Ja Jb- 



(3) 



Since oxide grows when metal ions reach the oxide/gas 
interface, the growth rate is 



dL 



(4) 



where U, is the volume per metal ion in the newly formed 
oxide. Comparison with ^ immediately shows that 



jb ^ t 



-1/2 



(5) 



for the asymptotic behavior of the ion fluxes. By the 
Planck-Nernst equations, the electric field E and the ion 
density gradients {dxTia, dxnb) also have the same behav- 
ior. Moreover, the quantities {ja,jb,E,dxna,dxnii) are 
all uniform throughout the oxide. 

We assume that the metal and the gas are neutral, so 
by Gauss's Law the surface charges (E(0), S(L)) and the 
electric fields {E{0),E{L)) are related by 



E{0)^ 



E(0) 



EiL) 



E(L) 



(6) 



Moreover, by continuity the assumption that there is 
charge and current only within the oxide leads to the 
conditions 



J(0) 



dE(0) 
dt 



J{L) = - 



dE(L) 
dt 



(7) 



This model immediately poses the questions of whether 
Mott's solution is self-consistent, and whether it is the 
beginning of an asymptotic series in powers of t^^/"^. 
To answer self-consistency, note that the uniform but 
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decreasing-with-time E ^ t~^/'^ leads to no bulk charge 
and to interfaces with equal and opposite charge, so they 
behave like a capacitor. Since E decreases with time, so 
must the charge on the capacitor. However, the model 
assumes zero current. Hence Mott's solution is not self- 
consistent. 

Nevertheless we will show, in response to the question 
about asymptopia, that each of the continuous variables 
can be expanded in an asymptotic series in where 
Mott's solution corresponds to n = 1, and that non-zero 
current J appears at order n = 3. We will also show that 
all of the continuous variables can depend upon position. 
This means that the bulk can develop a local and total 
charge density, with the surfaces not having equal and 
opposite charges, so that the capacitor model holds only 
to lowest order. It is already known that steady non- 
equilibrium current flow can cause local charge densities 
in the bulk.EH 

The fact that an expansion can be made of jb in powers 
of ^-1/2 igads to 



II. ON TWO-COMPONENT TRANSPORT 

Let the subscripts a and b denote electrons and metal 
ions respectively. Let no and uq/Z be the uniform equi- 
librium concentration of electrons and metal ions respec- 
tively, {na,nb) be their additional concentrations, {va,Vb) 
be their mobilities, {Da,Db) be their diffusion coefficients, 
(qa — — e, Qb = Ze) be their charges, E be the electric 
field, ks be the Boltzmann constant, T be the tempera- 
ture, and X be the position. 

Note that we use the Einstein Relations to rewrite the 
mobilities, which can have either sign, in terms of diffu- 
sion constants, which are always positive. Since we al- 
ways consider Z positive, if we want to consider oxygen 
ions and holes as the carriers, then only the sign of the 
electric charge e must be changed. For M^"*" and 0^~, we 
let e — > 2e and Z = 3/2. Therefore our results are quite 
general. 



A. Equations for Two-Component Transport 



so that 



Bint 



(8) 



(9) 



to higher accuracy than given by the pure parabolic law. 
This prediction can be subjected to experimental study. 
In practice, we must assume that 



L = V2At^/^ + B\nt + C + .. 



(10) 



because \nt is of order unity. We have not found data of 
a high enough precision to verify this In t correction. 

Section II gives the five fundamental equations for 
the five continuous variables (two continuity equations. 
Gauss's Law, and two Planck-Nernst equations), and in- 
dicates the expansion in powers of (The fact that it 
is such an expansion relates to the leading order term and 
the fact that the Planck-Nernst equation is a nonlinear 
function of the density.) Section HI gives the n — 1 so- 
lution, which has two interface-reaction-rate-determined 
integration constants: both A and a constant-in-space 
density not included in the Mott solution. Section IV 
presents the results of the n = 2 solution that gives the 
first correction to Mott's solution, and shows that there 
is a non-zero bulk charge density and a corresponding 
spatial variation to the field. There are four interface- 
reaction-rate-determined integration constants for n = 2. 
All higher-order solutions involve four interface-reaction- 
rate-determined integration constants. Section V consid- 
ers the oxide thickness growth rate. Section VI provides 
a summary and conclusions. The Appendices explicitly 
give the n = 2 and n = 3 solutions. 



We will employ the Einstein Relations 



Da " fesT' Db 



For electrons and metal ions, —Zqa = qt, so 

1 

Vt' 



VbDa = -ZvaDb, -Tr = i^ 



where 



Vt = 



(11) 



(12) 



(13) 



denotes a thermal voltage. Mott and Cabrera's papeif^ 
used {va,Vb) to denote electron and ion mobilities, which 
they may have intended to be (i^a, i^b)- 

The one-dimensional Planck-Nernst equations for the 
number flux densities associated with metal ions and elec- 
trons are 



ja = l^a{no + na)E - Dad^Ua, 

jb = Vb (J^ +nb^ E - Dbd^Ub. 



(14) 
(15) 



Rewriting mobilities in terms of diffusion constants using 

(fTTl), 



ja = -77^(flO + na)E - DadxUa, 

Vt 

jb = 77^ ("0 + Znb)E - Dbd^Ub, 

VT 



We also use the number continuity equations, 
dtria + dxja = 0, dtUb + d^jb = 0, 
and Gauss's Law, 



d^E 



-{Zub - ria), 



(16) 
(17) 

(18) 
(19) 
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where e is electron charge in Coulombs, and e is the per- 
mittivity of the oxide. 

With the charge density and current density defined 
by 

p ^ -e{na - Zub), J = -e{ja - Zjb), (20) 

use of the number continuity equations yields the charge 
continuity equation 



dtP + d^J — 0. 



B. Expansion Notation 



(21) 



We seek a series solution in powers of t^^^^ for the elec- 
tron and ion concentrations (densities) and fluxes, as well 
as for the electric field. They must satisfy the Planck- 
Nernst equations, the continuity equations, and Gauss's 
Law. Following Mott, we take the lowest order fluxes and 
field to vary as t~^/'^. Examining the structure of the 
Planck-Nernst equation for electrons (161, and inserting 
terms of order the nonlinear term UaE will contain 

terms of order t~^. Iteration yields that the series must 
be in powers of i^"/^ for integer n. 

We thus make an expansion of the form 



ja ^ ^ Jan^ ^7 jh — ^ ^ Jbn^ 



-n/2 



^a = ^ Nant" 



n/2 



rib 



n/2 



n=l 



-n/2 



-n/2 



-n/2 



(22) 
(23) 
(24) 
(25) 
(26) 



n=l 



n=l 



Here Jan, -hn, Nan, Nbn, En, Pn, and J„ are functions of 
the position along the direction of growth, x. From the 
above definitions, the dimensionality of {Jan, Jbn) is con- 
centration times velocity times i"/^, the dimensionality 
of {Nan, Nbn) IS Concentration times t"/^, the dimension- 
ality of the surface charge density (eI^-*, sl^') is charge 
per area times t"/^, the dimensionality of £'„ is electric 
field times t"/^, the dimensionality of p„ is charge den- 
sity times t"/^, and the dimensionality of J„ is current 
density times t"/^. 



C. On Specifying Chemical Reaction Rates at 
Surfaces 

In the presence of a true chemical reaction at a sur- 
face there is a single reaction rate, typically specified by 
a Butler- Volmer relationj^SEi] between the fluxes of all of 
the relevant components. In the present case the fluxes of 



the carriers are independent of one another, so that there 
are two statements about carrier fluxes at each surface, 
for a total of four conditions. With ft denoting an electro- 
chemical potential, near equilibrium (as we have here, in 
the asymptotic regime) each flux j will be proportional to 
its corresponding Aji across the interface (either M/MO 
or MO/0). The proportionality constant will depend de- 
tails of the reaction and the baseline properties of the 
system. That is. 



ja,b — Ga,bAp'aM 



(27) 



at each surface, so there are four G's. The equation 
becomes nonlinear far from equilibrium. Thus the ja,b 
are proportional to a non-equilibrium quantity, which we 
take to be a field Ei , as in Ref . [5l All of the unknown 
integration constants will be linear or higher in Ei. 

We will not attempt to carry this procedure any fur- 
ther. It is sufficient for our purposes to know that this 
can be done, and that in the present problem there are 
four constants associated with boundary conditions at 
the two surfaces for the two carriers. In principle, all 
of the quantities appearing in the solutions to the trans- 
port equations are determined by these surface reaction 
rates. For a true chemical reaction, which we expect to 
be described by a Butler- Volmer equation, the fluxes at 
each surface, because they are related, will be described 
by only a single independent coefficient G. Note also 
that the Butler- Volmer equation is non-linear, so that 
the boundary conditions can be nonlinear. Because we 
do not consider the boundary conditions in detail, we will 
neglect this possibility. 



III. RELATIONS BETWEEN EXPANSION 
COEFFICIENTS: ALL n 

A. Continuity Relations and Charge Conservation 

The continuity equations imply charge conservation, so 
we treat them in the same subsection. 



The continuity equations ( 18 1 yield 



n— 1 m— 1 

(28) 

for subscript a, and a similar relation holds for b. With 
m — (n — 2), so that J2m=i ~^ X]n=3' comparison of like 
powers of t yields, for n — 1 and n — 2, 



^xJan — 0; dxJbn — O5 

and, for n > 3, 



(n = l,2), (29) 



dr. J a 



n - 2 



2 



N, 



a(n— 2) , 



Nb{n-2), 



{n > 3). (30) 
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By definition we have 

Pn = e [ZNbn - Nan) , Jn= e (ZJ^n - Jan) , (31) 

and charge conservation for each n is 

dtPn + d^Jn = 0. (32) 
Charge conservation at each surface yields 



dt 



dt 



so 



i:(-i) 



~J\x=o = -eiZjb - ja)U=o, 



j.(0)^-(n+2)/2 



(33) 



= -e ^ (ZJb„ - Jan) |(.=o)t""/', (34) 

n=l 

n=l 

= e ^ (^>^()n - -^an) |(x=L)i""^^- (35) 
n=l 

Comparing powers of t we have, for n — \ and n = 2, 

^^f)ri|(3;=0) = Jan\{x=Q)-, 

ZJbn\(x=L) ^ Jan\{x=L), ("-=1,2), (36) 

and, for n > 3, 

~ ^ ^ ^ ^ ^'^{ZJbn — Jan) \{x=0): 

- ^ ^ ^ (^"^fen ^ "^an) |(x=L), 

(n>3). (37) 

Charge conservation over both surface and bulk yields 



n=l 



'/ y^iNan- ZNbn)t-''/^dx, 
■^0 n=l 



(38) 



so that, for each n, 



E(o) + Si^) = e / (iV™ - Zmn)dx. (39) 
Jo 



B. Gauss's Law 



Gauss's Law ( 19 ) reads 



n—1 \n—l 



Pnt 



-n/2 



(40) 



so, for each n, 



dxEn = (41) 
Gauss's Law at the surfaces ^ gives, for each n, 

y(0) y{L) 
En{0)^^, En{L)^ ^. (42) 



C. Planck-Nernst 



The Planck-Nernst equation for species a, ( 16 ), can be 
written as 



^ ^ J ant 



-n/2 _ _ 



Vt 



Uq 



Y^Ent' 



V2 



Da 

Vt 



J2 ^ani?™i- - Da Y,{dxNa)t~^/\ 

n=l 

(43) 



m.n— 1 



5]('^an + ^noEn + DadxNan)t-'''^ 



Vt 

Da 

'Vt 



an-^m^ 



)/2 



(44) 



m,n— 1 

with a similar form for species h. 

Db 
Vt 

ZDb 



J2(Jbn - ^noEn + DbdxNbn)t''''^ 

J2 iV6„i?„,t-('"+")/2. 



(45) 



Vt 



By matching coefficients of powers of t, we obtain the 
resultant equations for any n. 



D. Solving the Transport Equations 

There are five first-order differential equations for five 
continuous variables, so there are five integration con- 
stants at each order. For n > 3 the current density J„ 
at a; = L is known from the n — 2 value of surface charge 
density Therefore only four integration constants 

need be determined. These can be thought of as fixed by 
the "reaction rates" of each charge carrier at each of the 
interfaces, which will not be specified. 

The cases n = 1 and n = 2 are a bit simpler than 
n > 3. Nevertheless, for each n our strategy will be 
the same: (1) from the continuity equations find the ion 
fluxes Jan and Jb„; (2) from all of the five equations find 
an equation for p„ and solve it; (3) use this pn in Gauss's 
Law to find En] (4) find Nan and Nbn by substitution of 
Jam Jbn and En iuto thc Planck-Nernst equations. 



IV. SOLUTION FOR n = 1 

For n ~ 1, the calculations are simple, but illustrate 
what happens in higher orders. The continuity equations 
(29) give Jai and Jbi to be uniform, and charge conser- 



vation at each surface ( 36 1 gives 



Jai — ZJbi. 



The Planck-Nernst equations ( 44 1 and ( 45 1 yield 

Jai + ^noEi + Dad.Nai = 0, 

Jbi - jj^noEi + Dbd^Nbi = 0, 
Vt 



(46) 

(47) 
(48) 



where we take Ei as a uniform, experimentally deter- 
mined value. 

From the uniformity of Ei, 



d^Ei = 0, 
so that Gauss's Law yields 

Nal^ZNbl. 



(49) 



(50) 



Substitution of fluxes from ( 46 ) and concentrations 



from ( 50 ) into the Planck-Nernst equations gives 



Jai + ^noEi + Dad^Nai = 0, 



Jai Db 



noEi + Dbd^Nbi = 0. 



(51) 
(52) 



Thus, with (50), we have 



Jal^ZJbl=-{l + Z) 



DqDb npEi 
Db-D, Vt ' 



here, 



iV,i = ZNbi = 



H 



HuqEi 
Vt 



x + Mi; 



Da + ZDi 
Db - Dn 



(53) 
(54) 

(55) 



is dimcnsionless, and Mi is a constant of integration 
(units of concentration times s^^^) determined by the 
surface reaction rates (and thus linear in Ei). (Recall 
that (Jai, ^ai) are first order coefficients, which must be 
multiplied by in order to find the flux and num- 

ber densities (j, n).) Although constants of integration 
associated with reaction rates were discussed earlier, we 
consider them more explicitly here. 

Figure [T] illustrates the effect of changing the value of 
Ml. If All were very large, the high concentration of 
metal ions near the metal/oxide surface would oppose 
new ions from entering, whereas the high concentration 
of metal ions near the oxide/gas interface would encour- 
age more ions to be deposited on the oxide/gas surface. 



Metal 



cions 
-J sec 




Oxygen 



.1 1 ions 
J- J- sec 



Tions 
/ sec 



FIG. 1. The effect of the constant concentration Mi of metal 
ions. Here A^^i has too high a value of Mi. Only the par- 
ticular value of Ml in A'^^'' permits an equal rate of ions to 
enter and leave the oxide, here taken to be 7 ions per second. 



Eventually, the number of metal ions in the bulk would be 
insufficient to maintain the high rate of ions exiting the 
oxide, and would drop to some equilibrium value. Thus, 
Ml is determined by constraining the oxide to have no net 
ion-loading or ion-unloading in the bulk at order n = 1. 
Note that Ei, which is proportional to (the parabolic 
growth rate coefficient), is also related to the surface re- 
action rates. In general, net surface reaction rates involve 
a Butler- Volmer equation, but not far from equilibrium 
(as in the Mott solution) they can be linearized in the dif- 
ferences of various electrochemical potentials. This will 
ensure that there is no net surface reaction rate in the 
limit of equilibrium. 



From Gauss's Law at the surfaces (42) 



£i(0) = ^, Ei{L) 

€ 



.(L) 



Since Ei is uniform. 



.(0) 



.(L) 



eEi. 



(56) 



(57) 



For \va\ » Wb\ (or equivalently in this case, Da » 
Db), Mott and CabrergPfind for monovalent ions that 



Jai = -2D, 



dNai, 
dx 



(58) 



our results can be shown to be consistent with this. 



V. SOLUTION FOR n = 2 

We here summarize the n = 2 results. (For the ex- 
plicit solution, see Appendix |Aj) Recall that all coeffi- 
cients must be multiplied by t to find the physical vari- 
ables {j,n,E). With the constant M21 (in units s/m^) 
determined by surface reaction rates (and thus linear in 
El), the second order flux density coefficients are 



Ja2 = ZJb2 = -(1 + Z) 



f DaDb 



\ZDb + Da 



M21, 



(59) 
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and there is no net charge flux, 

•h = Ze.h2 - eJa2 = 0. (60) 

With the constants A/20, Pat^ and (each in units 
of s/m) determined by surface reaction rates (and thus 
Unear in Ei), the n = 2 coefficients of the concentrations 
of electrons and ions are given by 

Na2 = M20 + M21X + Pit^e^"' + -Pi2"^e-^/'= , (61) 



^62 = 



Mo, 



1 HtEl M21 



Z Z Vre 

^a2 * 

Here, Is is the screening length 

h = 



^a2 ^ ^a2 ^ 




(1 + Z)noe2 



(62) 



(63) 



There is a net charge in the bulk, given by p2t ^, where 



P2 



Vt 



(64) 



and 



Pi+) = -(1 + ^)ePi2+\ Pt'' = -(1 + Z)ePi-\ (65) 

Note that p2 has, in addition to surface charge within 
a screening length of the two surfaces, a uniform charge 
density with sign determined by e/{Da — Db) and inde- 
pendent of the sign of Ei (or, equivalently, the direction 
of current flow). Since Da ^ P'b here, the term is pos- 
itive. For holes and oxygen ions the carriers, we have 
—e/{Da — Db), but Db Da, so it is again positive. 
As found in previous work,^^ this uniform charge density 
leads to a quadratic voltage profile within the bulk, not 
within a screening length of either surface. 

The second order coefficient of the electric field is 



E' 



HEl 
Vt 



VtM2i MiEi 



noH no 
The surface charge coefficients are given by 



.(0) 



P^) - Pi- 



no 

eHEf 



(66) 



(67) 



Vt 



-L 



€VtM2i eMiEi 
noH no 

We have verified that 



s(°^ + + I p2dx = 0, 



(68) 



(69) 



so there is no net charge in the system. 



VI. RATE OF GROWTH OF OXIDE LAYER 

The oxide layer grows as metal ions reach the MO/0 
surface, and are taken into lattice positions to form a new 
oxide layer. Thus, the rate of growth of the oxide depends 
on the rate jb at which metal ions arrive, according to 
Q, or dL/dt = rijb. Including n = 3 (see Appendix [B ) 
we find the metal ion number flux (using Eqs. (53 1, (59) 
or ([ATt]), and 



Jb 



DaDb no El 



1 



Db 
Z 



Z 

fHnoEi 
V WtZ 



Da Vt 

DaDb 



ZDb 

X - 

2Z 



-Da 

Z 



-1/2 



M21 j t 

eEh\ 
2Ze J ' 



-3/2 



(70) 



where K3 is a constant of integration (units of flux den- 
sity times s'^/2) determined by interfacial reaction rates 
(and thus linear in Ei). Note that, if Db < Da (as for 
ions relative to electrons) then the microscopies must give 
i?i > for a positive growth rate. Keeping only terms 
of second order, integration of ( 70 ) with respect to time 



gives 



where 



A = 2 
B 



L = V2At^/^ + B\nt + 



l + Z^Y f DaDb \ nlEfn 



(71) 



Z 

DaDb 



Db-D„ 



V^ 



(72) 



D, + Z^Db 



(l-^)^-^M2r|a 



Vt 



Z 



(73) 



VII. SUMMARY AND CONCLUSION 



We have shown that the approach taken by Mott for 
parabolic growth of oxide films can be turned into a 
consistent asymptotic expansion, and we have explicitly 
given the form of the lowest three orders. Up to four in- 
tegration constants appear at each order, related to the 
surface reaction rates. At higher order the bulk film is 
found to be charged, with a corresponding non-uniform 
current density. 

The Appendices present the n — 2 and n = 3 solutions 
in detail, to show that the method can be used for any 
n to find the fluxes, concentrations, surface charges and 
electric field. As a consequence one can have confidence 
that the Mott solution gives the leading term in the com- 
plete solution of the complete set of transport equations. 

The most easily verifiable prediction from the view- 
point of experiment is the prediction that the first correc- 
tion to the linear growth law is logarithmic. Because Int 
is of order unity, data should be analyzed with an addi- 
tional constant: L = \/2At + B \nt + C + . ... A sampling 



7 



of the current literaturel^Hlll did not find enough preci- 
sion to confirm the logarithmic form of the correction 
term. 

We would like to thank Allan Jacobson for his com- 
ments and suggestions. This work was partially sup- 
ported by the Department of Energy through grant DE- 
FG02-06ER46278. 



Appendix A: Explicit Solution for n — 2 



By the continuity equations ( 29 1 Ja2 and Jb2 must be 



uniform, and by ( 36 ) they are related by 

Ja2 = ZJi,2- 



(Al) 



For n = 2, the Planck-Nernst equations from (44) and 



(45) are 



Ja2 + ^noE2 + Dad^Na2 = -^NalE^, (A2) 

Jb2 - ^noE2 + Di,d^Nb2 = ^NbiEi. (A3) 



Taking spatial derivatives of (A2) and (A3), and using 



Gauss's Law (41) and the uniformity of Ja2 and Jb2, 
yields 



- -f7 — P2 + Dbd^Nb2 = -r^d^{NbiEi] 



The right hand sides are found from ( 54 ) ; 

— P2 + Ua(Jx^^a2 — -^a i 

Kye Vrp 

- -y^P^ + DbdxNb2 = Db ^2 , 



(A4) 
(A5) 

(A6) 
(A7) 



where H is defined in ( 55 ) 



Subtracting (A6) multiplied by l/Da from (A7) mul- 
tiplied by Z/Db yields 

The solution to this equation, with 

y (1 + Zjnoe y (1 + Z)nQe^ 

and with new integration constants ^2^'' ^2 ^ 

r2 



P2 



HtEj 



(AlO) 



From (AlO) we infer that A*'a2 and ZNb2 are polynomials 



whose terms that are linear or higher are equal, and they 
may have different exponential terms. 



Substituting p2 from (AlO) into Gauss's Law (41 ), and 
integrating dxE2 yields 



Eo — 



Is 



R 



i+)pX/h _ pi-)^-x/l 



HEl 
Vt 



F2, 

(All) 



where F2 is a new integration constant with units V-s/m. 
By (54), {Nai,Nbi) are linear in x, so that the right- 



hand-side of the Planck-Nernst equation ( A2 ) is linear in 



X. Moreover, the continuity equation ( 29 ) implies that 



Ja2 is constant. Therefore by (All) for E2, the Planck- 



Nernst equation allows dxNa2 to be linear, so Na2 can 
be quadratic in x. Moreover, the exponential terms in 
{Da/VT)noE2 + DadxNa2 must cancel. Finally, from 
(AlO) any linear or quadratic terms in iV£i2 and ZNb2 



must be equal. With (M20, Af2i, M22, P^^'') being new 
integration constants, we therefore conclude that the fol- 
lowing form must hold: 



Na2 = M20 + M21X + M22X -r ^ „2 ^ ^ a2 



Use o f (l3lj) for n = 2 (p2 
(|A10| gives 



2^p( + )g.//. ^pMg-./i._ 

(A12) 

e{ZNb2 - Na2)), and p2 of 



Mo 



b2 



1 HeEf A/21 



M, 



Z Z Vtc 



Z 



P. 



(+) 



eP, 



(+)' 



a2 



PI 



z 

{-) 



22 2 
X 



eP, 



a2 



Ze 



Ze 



-X/Is 



(A13) 



Addition of ( A2 ) divided by Da and ( A3 ) divided by 
Db gives, with (Al) and ([54]), 



Ja2 



1 



1 



Da ZDb 



dx {Na2 



Nb2) = 0. 



(A14) 



Substitution for A*'a2 and Nb2 from (A12) and (A13) al- 



lows us to solve for some of the constants, 



M- 



21 



P 



1 

Z 

(+) 



2M: 
1 

^Z 



22 



P 



1 

z 

{+)' 



? p 



(-) 



1 + 



P 



Ze 



-x/h 



— —Ja2 



1 



Z^Dt 



(A15) 



So, since Ja2 is uniform, comparison of powers of x yields 
the conditions 



p(+) 
-'2 

Ja2 — 
M22 -- 



-{l + Z)eP^l 



(+) 



ZJb2 

-- 0. 



-(1 + ^) 



D„Db 



-(l + Z)ePi2'\ 



ZDb + Da 



M21, 



(A16) 
(A17) 
(A18) 
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We may thus rewrite (A12| and (A13) as 



a2 



Integration of these two equations gives two integration 
constants. However, these two constatnts are constrained 
M20 + M21X + P^2^e^/^= + p'j^'e"'^/'^ (A19) charge conservation across the interface at x = (37l, 



Nb2 = 



'J- 20 



1 HeEl 



Z Z Vre 



M21 



pi+)„x/U _ pi-)„-x/h 



a2 



(A20) 



Note that in ( A2 1 the exponential and linear coeffi- 



cients already match, by construction. A new constraint, 
however, is found by comparing the constant terms. 



DaDb 



Z^Di, + D 
Vt 



A/21 



DaUQ 

Vt 



DaM21 

(A21) 



(0) 



ZJb3\x=0 — '^a3|a;=0 



2e 



(B2) 

where Y.^ is given by (57 1. Integration of the continuity 
equations (Bl| then yields only a single new integration 
constant, which we call if 3: 



WtZ ^ 2z''^ Z +2Ze- 



(B3) 
(B4) 



There is a uniform net electric charge flux (i.e., current 
density) at order n = 3, 



yielding the condition that 

¥rM2i Ml El 



Fo = 



noH 



no 



We are thus left with four independent constants of in- 
tegration, which we take to be M20, M21 , P2^^ , and P2 ' . 
The reaction rates for electrons and ions at each inter- 
face, needed to produce the correct ion fluxes, provide 
the 4 conditions necessary to solve for these constants. 

For completeness we note that, from Gauss's Law at 



each surface ( 42 ) 



.(0) 



E2{0) = 



E2{L) = -- 



(A23) 



J3 = e {ZJb3 - Jqs) 



eEh 
2 ' 



(B5) 



(A22) due to the discharge of the surfaces in order n = 1. 



For n = 3, the Planck-Nernst equations from (44 1 and 



( 45 1 are 



- iNa2Ei + NalE2) , 

(B6) 



Jb3 - ^noEs + DbdxNb3 = ^ (iV^^^i + NbiE2) . 
Vt Vt 



(B7) 



Taking spatial derivatives of (B6) and (B7), and using 
Gauss's Law (41 1, 



Then (|All|) and (|A22|) yield 



.(0) 

^2 



L P 



.(+) 



eVrAhi eMiEi 



-L P 



P 



noH 



no 

eHEl 
Vt 



(A24) 



L 



eVtM 21 
noH 



eMiEi 
no 



(A25) 



d.Ja3+^P3+DadlN,3 



Vre 



D. 



,,-dANa2Ei+NaiE2), (B8) 
Vt 



dxJb3-^P3 + DbdlNi 



VTt 



b3 



ZDi 
Vt 



dANb2Ei + Nt,iE2) . (B9) 



Appendix B: Solution for 71 = 3 

Recall that in order n — 3, all coefficients must be 
multiplied by t'^/^ to find the physical variables {j, n, E). 

Using the n = 1 coefficients for electron and metal ion 
number densities from ( 54 ) and the continuity equations 
from (30) for n ~ 3, 



HnoEi Ml HnoEi Mi 

d^J^3^^^X+—, dxJb3 = ^^X+—X. 

(Bl) 



Subtracting ( B8 ) multiplied by 1/ Da from ( B9 1 multi- 



plied by Z/Db yields an equation for pa, 



e,--i-)a.J.3-(l + ^);|^P3 + a^P3 



1 
J 

eE 



TT^{dxNa2 + Z^dxNb2) + —{d^Nal + Z^dxNbl)E2 

Vt Vt 



Vt 



{Nal+Z^Nbl)d.,E2. 



(BIO) 



Substitution for (iVai,A'^bi) f rom ( [54| ), (A^a2,A^b2) from 
( |AT9l ) and (|A20|, E2 from ( |ATT| ), and J^g from (|B3| 
yields 
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f Db^-Da\ { HnoEh , Mi \ 
^[ D^D, )[ 2Vt 2 J 



The solution to this second order differential equation, From (B12) we infer that Nas and ZN^^ may include 
with two new integration constants (3['^^ and l3['\ is polynomials whose terms that are quadratic and higher 

are equal, and may include exponential terms that differ. 



P3 



(B12) 



where, with substitution of F2 from ( A22 1 



ai 



a2 = 



,(+) 



?(-) 



,(+) 



DaDb J 2 

Db-Da\ eHEi 



no 



HoVt 
(B13) 



DaDb J 2(1 + Z) 
1- Z^\ El HEi 



, Ml 



(B14) 



l + Z J 2Vt Wt 



2noL 
Ml 



,(+) 



(B15) 



1 - Z2 \ Si HEi 
l + Z ) 2Vt ~^ ~ 2nols 



>(-) 



g^i p(+) 



p(-) 



AVtIs 



(B16) 
(B17) 



Substituting from (B12) into Gauss's Law (41 ), and 
integrating d^E^ gives 



e 



Is 



"2 2 

2e 



.-X/Is 



2hpi+A 



- 2/,/3. 



(-) 



(B18) 



where F^ is a new integration constant, with units V- 
s3/Vm. 

Substitution of the coefficients Nai from ^5^, Na2 



from (A12), E2 from (All), J^s from (B3|, and E3, from 



(B18) into the n = 3 Planck-Nernst equation for elec- 



trons (B6) yields 



10 



M20E1 , M1F2 . noFg , K3 



Vt 



Vt 



Vt Da 



HMiEl M21E1 HnoF2Ei noai Mi 



Vt 



V^ 



eVT 2Da 



H^rioEf nQa2 HhqEi 
V^ 2eVT ^ 4:DaVT 

EiR 



(+) 



(1 + Z)eVT 



E1P2 



eVi 



(1 + Z)eVT eVr 

(+) 

~ eVT 



.Vt 



lJ2-^+2fA-^)]e-^/^^ 



eV^ 



X/Is 



eVT 



xe 



(B19) 



one new integration constant Afao, is 



+ 7^^:Ee--/'= +7^^^a;2e-/'= +7^^a;2e--/'= 



(B20) 



7^^ 



7^^ 



7^^ 



where 








M31 = 




(M20E1 


M1F2 


- 


\ Vt 


Vt 


M32 = 




( HMiEl 


M21E 




K V} 


Vt 



M33 = 



7^^ 



3 I V^ 2eVT 4DaVT ) ' 



{H + l)hEi Ml \ P 



.(+) 



Ft no J {1 + Z)e 

(-) 



(1 + Z)e 

Ml ^ {H + l)lsEi\ F2 
no 



Vt 



{1 + Z)c 



HEiP, 



(B21) 

(B22) 
(B23) 

(B24) 



{l + Z)e 

^ + 2?,/3(-) + Qfjf^) , (B25) 

(1 + Z)eVT (1 + (^^""^ " ^^'^^''O ' 

(B26) 



(-) 



HEiP^ 
'JlTz)eVT {l + Z)e 



(1 + Z)e' 



73 



(-) 



(1 + Z)e 



(B27) 
(B28) 



Use of 131]) for n = 3 {p^ = eiZNbs - iVas)), P3 from 



11 



(B12), and Nas from (B20l gives 





+ eA/30 


V 


Ze 


A/32 
7 


).^ + 
/ 




+ nt^ 




Ze 




+ 672^' 




Ze 






Ze 



02 + eM^i 
Ze 



Ze 
Ze / 



xe 



Ze 



(B29) 



We now have J ad,-, Jbs, E^, Na3 and Nii^ with five con- 
stants of integration, iir3, l3['^\ P\~ \ F3 , and M30. To 



reduce this from five to four, we add (B6) divided by Da 
and (|B7|) divided by D^, 



Ja3 Jb3 ^ ^ , Ar ^ ^ f /'2/?l , P1E2 
J^ + Tr+ dx(Na3 + Nb3) = — + 

Da Ub Vt \ e e 



(B30) 



As for n = 2, in (B30I the exponential, quadratic, and 
hnear terms match, by construction. A new constraint, 
however, is found by comparing the constant terms. 



K3 K3 eEi a2 + eM3i HeEf 

I TTT^ — I tt:^ — — r A131 + 



Da ZDb 2ZeDb 



Ze 



V^e 



(B31) 



Substituting M31 from (B21) yields 
Db 



F3 



D„ 



K3 



M20E1 M1F2 



{l + Z)DaDbJ no 

1 \ fZHeEf ^ eEi 
2D~b 



{1 + Z)ej V 



no 

Q!2 



(B32) 



a relation between integration constants F3 and K3. 

We are thus left with four independent constants of 
integration, which we take to be K3, M30, (3['^\ and /3[~^ . 
The reaction rates for electrons and ions at each interface, 
needed to produce the correct ion fluxes, provide the four 
conditions necessary to solve for these constants. 

For completeness, we note that from Gauss's Law at 
each surface (42), 



£^3(0) = 
Then, ^BT^ yields 

eF, 



E3{L) = -- 



V - -3 



(+) 



(B33) 



(B34) 



4^' = -.F3-aiL-^L^ 



(+) 



lsf3. 



- Is (M"^ 

+ ls (/3^'+2/,/?^"^)Le-^/'= 
We have verified that 



L/h 



.(0) 



psdx = 0, 



(B35) 



(B36) 



so there is no net charge in the system. 
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